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FORWARD 


The  eneloeed  ■sterlal  consists  of  lecture  notes  prepared  for  A  Short  Course 
on  Inertial  Ouidance  presented  by  Engineering  Extension  and  Physical  Sciences 
Extension,  UnlTeralty  of  California,  Los  Angeles,  October  10-21,  I960.  The 
theory  of  inertial  Instnnaents,  platforas,  and  error  analysis  are  presented 
in  other  parts  of  the  course.  This  oaterlal  is  concerned  nalnly  with  the  proper 
utilisation  of  Inertlally  derived  position  and  velocity  data  In  a  way  vhlch  vlU 
fulfill  the  aissloc  objectives. 
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ROCKIT  VKCLJ  OOISAffCE  SCHB<E8* 


1.0  DUMDucnoif 

The  guidance  scheae  or  philosophy  for  a  rocket  or  ballistic  Tehlcle  differs  fro« 
that  of  a  cruise  rshlele  In  tvs  principal  respects.  The  first  Is  a  result  of  the 
fact  that  a  rocket  Tehlcle  thrusts  at  high  lerels  for  relatively  avail  fractions 
of  Its  total  flight  tima,  ahlle  essentially  the  opposite  Is  true  for  a  cruise 
vehicle.  Since  significant  control  forces  are  available  only  during  povered  flight, 
the  rocket  vehicle's  guidance  syatsv  mat  be  able  to  direct  the  vehicles  course 
during  this  tlae  In  a  way  which  will  precisely  Influence  Its  position  and/or 
velocity  vlnutes,  hours,  or  even  vonths  later. 

The  second  major  difference  Is  In  the  area  of  diaenalons.  Rocket  guidance  systems 
are  three  dlvensioaal,  with  vehicles  being  guided  between  two  points  la  Inertial 
space.  Cruise  vehicle  guidance  systems,  on  the  other  hand,  are  basically  two 
dimensional  la  nature,  with  the  third  dimension  being  fixed  or  externally  supplied. 
Vehicles  are  guided  over  a  known  surface  like,  for  evaaple,  the  earth. 

1.1  Vehicle  Characteristics 

Rocket  vehicles  have  lengths  which  rsnge  froai  a  few  feet  In  the  ease  of  short 
range  missiles  to  100  feet  or  more  In  the  case  of  ICM's  and  space  vehicles. 

Weights  again  range  fron  a  few  pounds  to  hundreds  of  thousands  of  pounds.  They 
can  have  any  number  of  stages,  one  to  four  being  the  most  eosaxmi  and  nay  be 
either  liquid  or  solid  propelled.  Thrust  accelerations  during  each  stage  usually 
vary  from  1.5*3  g's  at  Ipiltlon  to  8-10  g's  at  burnout.  The  structu.*,  particularly 
for  the  larger  vehicles.  Is  extremely  ll^t  with  at  least  90  percent  of  the  gross 
welf^t  being  fuel  and  oxldlser.  The  resulting  lack  of  rigidity  often  places 
rather  severe  llmltatloos  on  the  vehicle's  maneuvering  capability. 

•The  author  would  like  to  ac^wledge  the  contributions  by  his  many  colleagues  in 
the  guidance  areas  at  Space'  Technology  Laboratories.  Re  would  la  particular  like 
to  single  out  J.  K.  Bachar,  P.  Baskin,  J.  A.  Joseph,  T.  U.  Layton,  J.  W.  McCarthy, 

V.  J.  McLaughlin,  R.  R.  Southworth,  and  D.  V.  Whltconbe  of  the  Inertial  Ouldance 
Department.  R.  M.  Page  of  the  Ouldance  and  Navli|atlon  Department,  and  V.  Schroeder 
of  the  Coaputers  and  Ouldance  Department. 


rrrL/Tn-6o-ocxx)-ORl86 

Page  2 


Flight  path  control  is  usually  obtalnsd  by  polntiixc  the  missile  Im  the  dlrec* 
tlon  of  desired  thrust,  vith  changes  In  attitude  being  obtained  by  momentarily 
deflecting  the  direction  of  thrust.  This  can  be  accomplished  by  gimbaillng 
the  thrust  chamber  in  the  case  of  a  liquid  or  the  nozzle  in  the  case  of  a 
solid.  Alternatively,  the  flame  pattern  can  be  deflected  by  means  of  Jet  vanes 
or  Jetavators  placed  in  the  thrust  stream.  Guidance  can  be  radio,  inertial 
or  a  coaibination  of  the  two.  The  diacussion  in  this  chapter,  vhlle  aimed 
primarily  at  inertial  guidance,  is  frequently  general  enough  to  apply  equally  veil 
to  other  types  of  systems. 

1.2  Trajectory  Characteristics 

\ 

The  trajectory  for  a  rochet  vehicle  can  be  divided  into  three  types  of  phases 
as  shown  in  Figure  1;  I  powered  flight,  II  free  flight,  and  III  re-entry. 


(a)  ICIK 


I  -  Powered  Flight 

Flight  Phases  \  II  -  Free  or  Balliatlc  Flight 


^III  -  Re-entry 


Figure  1  -  SCME  TYPICAL  ROCKET  TRAJECTC«IES 


5TL/T!«-60-0000-<Ha.86 

Page  3 


There  will  be  one  or  more  each  of  the  powered  and  free  flight  phases.  There 
will  be  a  re-entry  phase  only  if  the  payload  is  retviming  to  earth  or  arriv¬ 
ing  nt  some  other  planet  having  an  atmosphere. 

The  initial  powered  phase  is  the  most  complex,  because  of  the  exit  stscsphere. 
The  trajectory  usually  begins  with  the  missile  rising  vertically  for  a  few 
seconds.  During  this  time  It  rolls  to  the  proper  heading.  The  vehicle  then 
executes  Its  pitch  moneuTer;  after  a  short  transient,  usually  called  transi¬ 
tion  turn,  a  gravity  or  xero  lift  turn  begins  and  continues  until  the  Bissile 
has  effectively  left  the  atawaphere.  The  gravity  turn,  which  is  accoaqjllshed 
by  causing  the  sdssile  to  thrust  always  along  its  velocity  vector,  ainimizes 
drag  effects  and  aerodynaoic  heating.  The  gravity  turn  is  usually  continued 
to  some  staging  point,  although  this  is  not  always  the  case, particularly  when 
there  is  only  a  single  stage.  Mter  leaving  the  atmosphere,  structural  con¬ 
straints  can  be  relaxed  and  a  more  arbitrary  attitude  profile  can  be  pre¬ 
scribed.  A  very  high  acceleration  vehicle,  however,  can  achieve  the  desired 
velocity  before  it  ever  leaves  the  atmosphere.  This  can  cause  significant 
steering  problems. 

Whan  thrust  has  been  terminated,  the  vehicle  begins  its  free  flight,  where 
gravity  is  the  only  acting  force.  The  free  flight  trajectory  lies  completely 
within  a  plane  which  contains  the  center  of  the  earth  and  .will  be  in  the  shape 
of  a  conic  -  either  an  ellipse,  a  parabola,  or  a  hyperbola,  depending  on 
whether  the  velocity  is  below  or  above  escape  velocity,  the  parabola  being  the 
limiting  case.  In  the  case  of  a  ballistic  missile,  the  ellipse  intersects  the 
earth  at  the  target.  Actually  the  earth's  oblateness  causes  the  trajectory  to 
be  non-planar  and  to  differ  slightly  from  a  true  ellipse.  Similarly,  the  in¬ 
fluence  of  other  celestial  bodies  on  earth  satellites  and  space  probes  keeps 
then  from  being  pure  conics.  The  thrust-coast  sequence  can  bo  repeated  essen¬ 
tially  any  number  of  times  depending  only  on  the  mission. 
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Any  rocket  vhlch  returns  to  earth,  such  as  a  ballistic  alsslle  or  a  maimed 
space  vehicle,  nu*t  finally  undergo  a  re-entry  phase.  Here  non-noaiinal  re¬ 
entry  conditions  such  as  winds  or  density  variations  can  also  contribute  to 
syste*  inaccuracy,  since  these  effects  are  usually  not  predictable  during  the 
boost  guidance  period.  If  re-entry  guidance  is  used,  then  these  effects  are 
eliminated. 

1.3  The  Guidance  Problem 


Hie  thrust  of  a  rocket  engine  is  a  cootplex  function  of  the  engine  (and  propellant) 
parameters,  the  air  pressure,  the  temperature,  the  vehicle  acceleration,  and  to  a 
lesser  extent,  a  variety  of  other  quantities.  For  a  given  set  of  engine  and 
vehicle  parameters,  any  desired  trajectory  can  be  synthesised  by  using  simulation 
techniques  on  a  high  speed  digital  computer,  providing,  of  course,  that  the  per¬ 
formance  limitations  of  the  vehicle  are  not  exceeded  and  provided  the  relative 
positions  of  vehicle  and  target  are  known.  On  the  computer  the  desired  trajectory 
Is  achieved  by  time  programming  vehicle  attitude  and  terminating!  thrust  at  the 
appropriate  time  or  times.  How  then  la  the  powered  flight  of  an  actual  vehicle 
controlled  so  that  it  too  accomplishes  its  desired  mission? 

A  Simple  Example--  A  simple  example  for  purposes  of  illustration  the  German 
V-2  missile,  developed  near  the  end  of  World  War  11*’.  This  mlsBlle  had  a  take-off 
weight  of  28,600  pounds,  a  thrust  of  59,600  pounds,  and  a  maximum  thrust  accelera¬ 
tion  of  6.U  g's.  A  typlcaLL  trajectory  is  shown  in  Figure  2. 


Figure  2  -  V-2  THAJECTOFT 

•See  for  example ,  ne’rt i ol  Guidance  for  Rocket -Propelled  Missiles",  by 

W.  T.  Russell,  Jet  Propulsion,  January  1956- 
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For  thi»  ctse,  the  Kinalle  had  a  range  of  230  miles  and  nominally  biimed  out  at 
70  seconde,  about  15  miles  downraage  vith  a  speed  (V)  of  about  6000  fps  and  a 
burnout  angle  (r)  measured  .Trom  the  horl  rental  of  about  5**°- 

A  Biatpie  scheme  for  guiding  such  a  vehicle  la  ogalc  to  program  attitude  and 
bximout  aa  functions  of  time  aa  is  done  in  the  simulation.  A  pitch  attitude 
programmer  is  installed  In  the  missile  and  its  output  compared  with  the  gimbal 
angles  of  an  inertial  platform,  the  difference  being  used  aa  the  control  system 
error  signal.  Yav  and  roll  can  also  be  programed,  probably  to  tero.  A  deck 
is  used  to  shut  off  the  engine  at  70  seconds  or,  alternatively,  only  enough  fuel 
Is  placed  In  the  missile  so  that  the  engine  bums  out  at  70  seconds. 

The  acbeane  Just  outlined  performs  very  well  as  long  as  both  the  missile  and  its 
environment  are  nominal.  The  impact  accuracy  is  limited  only  by  the  performance 
of  the  programer  and  the  attitude  reference.  A  nominal  missile  is,  however, 
ocly  the  average  of  an  ensetnble,  vith  any  given  vehicle  differing  from  the 
nominal  to  some  degree.  Some  of  the  more  significant  pertuxhatlona  aa  far  aa 
the  trajectory  Is  concerned  ai’e  given  in  Table  1. 


Table  1 

Significant  Perturbations 

Thi*u8t  Variation 
Initial'  Maas” Varlatl  on 
Mass  Plov  Rate  Variation 
Thrust  Misalignments 
Drag  Variations 
Wind  (gusts  and  shear) 


Thrust  and  mass  variations  can  be  in  the  neighborhood  of  a  few  percent,  mis¬ 
alignments  around  a  degree,  and  winds  sometimes  in  the  hundreds  of  nrph.  While 
It  Is  clear  that  such  perturbations  can  cause  impact  errors,  it  Is  necessary 
to  examine  the  equations  of  motion  if  a  more  quantitative  indication  of  accuracy 
Is  desired. 
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An  Error  AnaJyolB  -  For  a  rnngc  of  230  miles,  the  earth  can  be  assumed 
flat  to  a  first  approximation.  If  x  la  dovnrange  and  t  Is  up  and  If 
the  re-entry  atasjsphere  Is  Ignored,  the  eqiiatlons  vhlch  describe  the  missile's 
flight  from  burnout  (BO)  to  Impact  are 
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where  x  ,  t  t  k  »  and  z  are  BO  positions  and  velocities,  t..  Is  the 

0  0  0  O  2 

time  of  free  flight,  and  g  Is  gravitational  acceleration  equal  to  32.2  fps 
The  problem  can  bo  further  simplified  If  It  Is  assumed  that  t  •  and  that 
x^  is  small  compared  to  x  Renee 
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From  (3)  and  (4)  It  follows  that 
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Prom  (6)  t.he  velocity  miss  coefficients  In  Table  2  are  easily  derived.  These 
coef flclontk''  relate  burnout  velocity  errors  to  Impact  position  errors. 
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The  effect  of  a  typical  perturbation  vlll  nov  be  exaadned.  Cooeider  the  caae 
where  the  thrust  la  51^  high  during  the  entire  70  seconds  of  powered  flight  ao 
shown  by  the  upper  path  In  Figure  3-  The  velocity  at  burnout  is  actually  the 


Thousands 
of  feet 
(Vertical) 


6000  fps 
5U° 


where  the  "T"  subseript  Indicates  thrust  and  the  "g”  subscript  gravity. 
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For  a  70  second  povored  flight,  i  is  efiuil  to  2254  fps,  assuadng  g  to 
¥e  m  constant.  lenee 

-  -05(3500)  -  175  fp«  (3 

-  .05(4800  +  2254)  -  353  fp*  (3 

Using  the  slss  coefficients  from  Table  2,  £x.  can  be  conputed  to  be 

dK  -  .•^175)  +  -044(353)  «  26.0  miles  (] 

Clearly  a  dassiraiige  miss  of  this  site  would  make  the  system  unacceptable  for 
most  appllcatloms. 


A  More  So>id.gtieated  System  -  If  the  system  Is  to  perform  adequately  In 

the  face  of  expected  perturbatloxu,  then  It  appears  that  some  closed  loop  scheme 

CBist  be  smplcyed.  The  Oeroans  took  the  first  step  In  this  direction  by  mounting 

an  Intsfrmtlng  aceeleroneter  along  the  V-2's  roll  axis.  Thrust  was  then  terminated 

ehen  a  pre-set  value  of  thrust  velocity  was  reached.  Again  consldef  the  case  of 

a  missile  with  thrust  51^  high.  Referrtng  to  (0)  and  (9)  it  Is  seen  that  the  only 

BO  velocity  perturbation  will  be  due  to  i  which  will  be  low  in  magnitude 

go 

since  the  tine  to  BO  Is  reduced  by  this  amount  and  g  has  been  assumed  con¬ 
stant.  The  reduction  In  time,  of  course,  la  due  to  the  fact  that  the  pre-set 
value  of  will  be  reached  earlier  due  to  the  high  thrust.  The  downronge 

miss  Is  now  costputed  to  be 


.044  dz 

o 

.044  (.05)(2254) 


5.0  miles 


Comparing  (13)  and  (12)  shows  that  the  downrange  miss  has  been  reduced  by  a 
factor  of  five.  A  further  reduction  could  probably  be  achieved  by  cutting  off  as 
a  function  of  both  time  and  thrust  velocity.  Beyond  this,  path  variations  In 
gravity  and  non-standard  burnout  position  would  have  to  be  considered.  Also 
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attention  vould  have  to  be  paid  to  the  lateral  direction  where  winds  and  mls- 
oilgnaents  can  cause  substantial  errors.  In  Section  2,  a  sore  basic  approach 
to  the  guidance  problem  is  taken  i  and  the  basic  elements  of  a  system  are  dis¬ 
cussed. 
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2.0  FliXENTS  OF  GUIDANCE 

In  the  previous  section,  It  was  shovn  that  a  considerable  improvement  in  system 
accuracy  can  be  obtained  with  a  relatively  simple  guidance  aystera  employing  a 
single  body  mounted  accelerometer.  It  is  easily  shown,  however,  that  the  future 
trajectory  of  a  vehicle  does  not  depend  on  any  single  variable,  but  on  all  six 
components  of  its  present  position  and  velocity  ( and  time  eilso  in  the  case  of  a 
moving  target).  It  is  therefore  obvious  that  in  general  any  really  precise 
g\ii dance  system  must  measure  and  utilise  (at  least  implicitly)  all  seve*.  vari¬ 
ables.  For  purposes  of  analysis,  the  system  can  be  divided  into  three  functions 

(1)  Navigation 

(2)  Computation 

( 3 )  Control 

In  the  paragraphs  that  follow,  each  of  these  functions,  as  well  as  their  combined 
capability,  is  examined  in  detail.  Inertial  components  are  assumed  to  be  perfect, 
since  component  error  effects  can  be  analysed  separately. 

* 

2.1  The  Navigation  Function 

Navigation,  as  defined  here,  consists  of  determining  the  vehicles  position  and 
velocity  in  some  known  frame  of  reference.  Related  to  this  is  the  computational 
coordinate  system,  which  can  be  either  inertial  or  earth  fixed.  It  can  be  of 
the  local  vertical  type,  but  usually  is  not.  Two  of  the  more  common  coordinate 
systems  are 

(1)  Launcb-Cetrtered  Inertial  -  The  system  is  inertial  and  is  centered  at 
the  launch  site  at  the  instant  of  launch.  It  typically  has  x  hori¬ 
zontal  and  in  the  launch  direction,  t  vertical,  and  y  completing 
the  right  handed  set.  It  may  be  desirable  from  a  computational  stand¬ 
point  to  rotate  x  and  s  somewhat  in  the  x-t  plane. 

(2)  Launch  Centered  Earth  Fixed  -  This  is  an  earth  fixed  coordinate  system, 
having  the  same  original  orientation  as  (l).  It  is  used  primarily  when 
it  is  desirable,  for  hardware  reasons,  not  to  remove  earth  rate  torqulng 
from  the  gyros  at  launch.  Computationally  it  has  both  advantages  and 
disadvantages  to  be  discussed  later. 

*It  is  possible  to  in  effect  combine  the  navigation  and  computation  functions 
into  a  single  operation.  Such  schemes  are  not  covered  here. 


STL/TH-o0-0000-(2ll86 
Pa«e  U 


Any  actual  Inprtlal  guidance  system  has  up  to  three  other  coordinate  systems 
of  interest.  All  four  have  the  same  center,  but  may  have  different  orienta¬ 
tions.  First  there  la  the  platform  coordinate  system  determined  by  the 
leveling  and  alignment  references.  Secondly  there  is  the  gyro  coordinate 
system.  Finally  there  is  the  accelerometer  coordinate  system.  OriRlnally 
there  was  only  a  single  system,  but  schemes  presently  being  used  to  reduce 
component  errors  and  simplify  computations  frequently  require  all  fo’jr. 

From  a  computational  standi>oint,  the  only  ones  of  Interest  are  the  accelero¬ 
meter  and  computational  coordinate  systems,  since  it  is  necessary  that  the 
airborne  computer  mechanise  the  api'roprlate  transformation  matrix. 

The  Itavlgatlon  Loop  -  Since  accelerometers  measure  only  non-gravltational 
forces,  the  total  acceleration  la  given  by 

•• 

R(t)  -  g(R)  J^(t)  (U)* 

where  R  la  inertial  position,  is  thrust  or  measured  acceleration, 

and  g  is  gravitational  acceleration.  The  block  diagram  for  a  navigational 
or  kinematic  loop  is  shown  in  Figure  U.  As  can  be  seen,  sensed  acceleration 
is  rotated  into  the  computational  coordinate  system  and  is  then  added  to  the 


Figure  k  -  THE  HAYIOATTCW  LOOP 


♦A  bar  over  a  variable,  as  in  g  ,  indicates  a  vector  quantity. 
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gravitational  acceleration  which  has  been  computed  based  on  the  cunputed  position. 
Position,  of  course,  it  obtained  by  doubly  integratiiig  the  total  acceleration. 

All  but  the  actual  sensing  of  acceleration  is  usually  accomplished  in  the  computer, 
in  most  cases  a  digital  contputer.  If,  as  is  frequently  the  case,  the  accelerometer 
actually  generates  pulses  representing  velocity  increments,  the  computer  accumu¬ 
lates  the  Increments  continuously  and  periodically  adds  the  total  to  the  integrated 
gravitational  acceleration.  When  accuracy  permits  an  analog  computer  can  be  used. 
Here  the  separation  between  sensors  and  computer  may  be  rather  obscure. 

The  presence  ofJUie, gravity  computation  is  one  of  the  basic  features  of  a  precise 
inertial  guidance  system  which  distinguishes  it  from  the  more  rudimentary  variety 
like  the  one  described  in  Section  1.  If  all  trajectories  were  close  to  nominal, 
there  would  actually  be  no  need  for  a  gravity  computation,  sines  the  effect  of 
gravity  could  be  pre-calculated.  Non-standard  missiles,  launch  delays  with  a 
non-earth-flxed  target,  winds,  etc.,  wi 11, however^ cause  non-ctandard  gravitational 
accelerations.  This  is  particularly  true  of  long  range  miBslles  where  both  the 
magnitude  and  direction  of  g  can  change  appreciably  during  the  long  powered 
flight  period. 

"Hie  Gravity  Computation  -  The  basic  expression  for  gravity  assuming  a 
round  earth  is 

8  “  -  R  (15) 

|Rr 

where  R  equals  the  position  vector  measured  from  the  renter  of  the  earth,  0 
is  the  universal  gravitational  constant,  and  M  is  the  mass  of  the  earth.  If 
oblateness  is  considered,  as  it  must  be  in  most  cases  ,  the  formula  contains  euldi- 
tlonal  terms.  From  a  mechanlratlon  standpoint,  there  are  two  objections  to  equa¬ 
tion  (15)-  First  the  coordinate  system  is  unnatural  for  an  inertial  system; 
second,  and  most  important,  the  expression  contains  a  division  and  a  square  root, 
both  of  which  are  relatively  slow  processes  in  airborne  digital  computers. 

The  above  problem  can  often  be  circumvented  by  observing  that  in  many  eases  the 
gravity  expression  need  be  valid  over  only  a  few  miles  in  y  and  a  few 
hundred  miles  in  x  and  s  If  such  Is  the  case,  (15)  can  be\expanded  in 

♦If  not  considered,  oblateness  cpji  cause  a  miss  In  the  order  of  10  n.mi.  for 
a  5500  n.mi.  ICBM. 
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a  Taylor  •erica  about  a  point  midvny  along  the  powered  flight  trajectory,  an 
example  being 

g  -  (fi  R^)  GM  (C^  +  C^x  >  C^t  +  C^x  C^^  +  C^xx)  (l6) 

Like  any  other  power  series,  the  number  of  ter**  required  depends  on  the  Taria- 
tions  in  and  the  accuracy  desired.  It  can  be  shown  that  all  of  the  coeffi¬ 
cients  are  functions  of  the  launch  latitude,  launch  arlauth,  and  expansion  point. 
It  is  possible  also  to  include  the  effect  of  centripetal  acceleration  in  the  ex¬ 
pression  when  an  earth  fixed  coordinate  system  is  used. 


Dynasde  Beha\d.or  -  Sotae  insight  into  the  dynamic  behavior  of  the  naviga¬ 
tion  loop  can  be  obtained  by  writing  the  components  of  (15)  in  terms  of  a 
rectangular  coordinate  system  with  its  origin  In  the  vicinity  of  the  trajectory 
and  its  X  axis  vertical.  Assuming  motion  in  the  trajectory  plane,  only  g^ 
and  g  will  be  considered.  These  can  be  written  \ 


-  GM  X 

-  aM(x  +  R^) 

~r^  71 

[x  (R^  +  z)  r' 

2 

where  it  is  easily  shown  that  GM  is  equal  to 
(14)  into  component  form  gives 


(17) 


(18) 


Expanding  equation 


Sc 


“Tx^  «x 


(19) 


(20) 


^6ee  for  exaaple,  "A  Polyncsiial  Gravity  Expansion  for  the  Anna  Inertial  Guidance 
System",  (Confidential),  by  D.  W.  Vfhltcoobe,  Space  Technology  Laboratories  Report 
«  41.2-37,  17  December  1956. 
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Writing  the  perturbation  equations  for  (19)  and  (20)  and  assuming  that  gravity 
perturbations  are  due, to  position  perturbations  gives 


da.  « 


sr  ^  * 

ST- 

'Si  *  sr  *'  * 

sr 

Taking  the  partial  derivatives  of  (17)  and  (l8)  In  the  vicinity  of  the  origin 
and  substituting  them  Into  (21)  and  (22)  gives  the  perturbation  equations  for 
the  navigation  loop 


Ax 

”o 


CH  -  —  ^  (24) 

The  solutions  to  (23)  and  (24)  for  constant  values  of  As^  perturbations  are 


n  \\ 


Az 


1) 


(26) 


Wear  the  surface  of  the  earth,  the  sinusoidal  oscillation  in  x  has  the  familiar 
Schuler  frequency  where  2s  »  64  minutes.  The  z  channel  i;;.  of  course, 

divergent  with  time.  For  the  three  dimensional  case,  the  behavior  of  the  y 
channel  is  Identical  to  that  of  the  x  channel.  The  point  to  be  br'^ught  out  Is 
that  the  computational  loop  is  "Schuler  tuned"  even  though  the  platf  rm  Is  not 
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locally  leveled.  Since ,  hovever,  rocket  vehicles  invariably  have  povered  flight 
tlsies  much  shorter  than  6U  minutes «  acceleration  errors  propagate  intc  velocity 
and  position  errors  proportional  to  time  and  time  squared;  from  a  perturbation 
standpoint,  the  gravity  loop  appears  to  be  open  for  time  Intervals  small  compared 
to  the  Schiller  period. 

Multiple  Quidanee  Periods  -  When  inertial  guidance  is  to  be  used  for  the 
launching  of  high  altitude  earth  satellites  or  for  space  travel,  additional 
povered  flight  periods  can  occur  after  long  periods  of  free  flight.  In  order 
to  establish  the  initial  conditions  for  the  navigation  loop  at  the  beginning  of 
a  subsequent  powered  flight  period,  it  is  of  course,  possible  to  operate  the 
navigation  loop  during  the  free  flight  period.  If  this  is  done,  (1)  accelero¬ 
meters  must  be  disabled,  since  any  output  during  the  free  flight  period  is 
erroneous  and  should  be  ignored,  (2)  the  gravity  formula  or  expansion  must  be 
valid  over  the  entire  free  flight  region,  and  (3)  the  digital  computer  may  con¬ 
sume  considerable  power^oyer .^-long  coast. 

A  second  method  for  establishing  the  initial  conditions  for  the  second  povered 
flight  or  bum  is  to  predict  them  from  the  observed  first  burnout  position  and 
velocity.  Since  the  guidance  equations  can  control  four  of  the  initial  condl- 
tlons,  it  is  necessary  to  predict  only  the  remaining  three.  The  prediction 
approach,  which  can  be  accoBgJllshed  with  either  explicit  formulas  or  polynomial 
expansions,  has  the  disadvantage  of  requiring  more  computer  meaory.  All  things 
considered,  though,  it  is  usually  the  preferred  approach  for  free  flints  of 
any  appreciable  duration. 

2.2  The  Computation  function 

Given  the  position  and  velocity  of  the  vehicle  as  well  as  the  tliae,  the  computa¬ 
tion  function  consists  of  utilizing  this  information  to  generate  error  or  status 
signals  ^ich  can  be  used  to  control  the  flight  path  of  the  rocket  in  a  way  which 
will  result  in  the  sdssions  being  accomplished.  The  functions  used  to  generate 
the  control  error  (or  status)  signals  are  commonly  called  "guidance  equations". 


♦T^is  is  shown  in  Section  2.2. 
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Thrust  Ma^itude  Control  -  Pocket  engines  whose  thrust  can  bo  '.ontrolled  in 
ma^ltude  as  well  as  direction  are  said  to  have  thrust  magnitude  control.  Vehicles 
with  this  feature  have  more  flexibility.  For  example,  they  can  bo  nade  to  fly  the 
exact  nominal  trajectory  and  hence  can  be  made  to  burnout  at  a  pre-specified  posi¬ 
tion,  velocity,  and  time.  The  guidance  computations  in  this  cane  cm  be  greatly 
•implified,  since  it  la  only  necessary  to  measure  the  three  cct.?onent8  of  thrust 
acceleration  and  to  compare  them  or  their  time  derivative  with  the  nominal,  profiles 
which  have  been  stored  as  functions  in  the  airborne  cosiputer.  By  controlling  the 
direction  (steering)  and  the  magnitude  (throttling)  of  the  thrust  it  is  possible 
to  match  the  stored  profiles  to  an  arbitrary  degree,  depending  only  on  the  response 
of  the  control  system. 

A  typical  conflgxiratlon  for  a  system  with  thrust  magnitude  control  la  shown  in 
figure  5-  Besides  being  simple  to  jsechanlae,  such  a  system  can  be  advantageous 


g 


Figure-^  -  A  OVIDANCE  STSTO!  WITH  THRUST  CC3WTR0L 


when  it  Is  necessary  to  match  a  desired  trajectory  exactly,  as  in  the  case  with 
the  rendeivous  problem.  It  should  be  pointed  out,  however,  that  the  combination 
of  launch  time  variations  and  a  non-earth-fixed  target  causes  the  nominal  trajectory 
concept  to  loose  most  of  its  meaning.  Throttling  in  any  case  has  not  been  con¬ 
sidered  desirable  from  a  propulsion  hardware  standpoint,  so  very  few  rocket  engines 
have  this  capability.  The  remainder  of  this  chapter  is  therefore  concerned  only 
with  the  class  of  vehicles  whose  thrust  can  be  controlled  In  direction  only,  and 
whose  thrust  magnitude  varlen  in  the  order  of  a  few  percent  from  engine  to  engine. 
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Available  Techniques  -  Most  of  the  available  techniques  for  guiding  rocket 
yehicles,  vere  developed  as  a  part  of  the  various  ballistic  missile  programs- 
Since  these  vehicles  are  intended  to  deliver  their  payloads  to  a  precise  point 
on  the  surface  of  the  earth,  the  first  approach  to  their  guidance  vas  to  control 
the  missiles  flight  p  .th  on  the  basis  of  predicted  errors  in  impact.  Functions 
of  the  fora 

I 

.  f^(R,  g,  t)  (27) 

-  f^(R,  R,  t)  (28) 

can  be  generated,  where  M.  is  the  downrange  miss  and  M  is  the  crossrangc 

a.  c 

mins,  both  at  the  target  altitude.  The  functions  will  differ  slightly  vlth  the 
miss  coordinate  system  employed.  The  one  most  often  used  is  the  so  called 
instantaneous  ii^act  pjint  (IIP)  coordinate  system.  Ilils  is  an  ox^hogonal 
coordinate  system  centered  at  the  target,  with  the  axis  defined  by  the 
downrange  IIP  locus  or  the  locus  of  Impact  points  which  occurs  when  only  the 
thrust  cutoff  time  of  a  nominal  vehicle  is  varied. 

Impact  error  in  the  actual  system  is  controlled  by  changing  the  flight  path 
(steering)  in  the  lateral  direction  (yaw)  until  M  is  driven  to  lero  and 
then  terminating  thiust  when  reaches  tero.  For  this  case,  pitch  steering 
is  not  required,  since  only  two  degrees  of  control  are  necessary.  Here  an  open 
loop  pitch  attitude  program  can  be  used,  e.g.,  const int  attitude. 

While  the  above  scheme  is  workable,  it  has  several  disadvantages.  First  of  all 

the  M.  and  M  expressions  are  both  functions  of  seven  vEurlables  and  are 
d  c 

laborious  to  compute.  Probably  even  more  Iwoortant  is  the  fact  that  the  expressions 
are  not  general  and  are  awkward  for  anything  but  ballistic  missile  (ICBM,  IRSM,  etc.) 
applications.  Renee,  the  "required  velocity"  concept  is  now  more  generally  used*. 


required  velocity  is  extracted  from, 


"An  Introduction  to 


Inertial  Ouldance  Concepts  for  Ballistic  Missiles",  (Tutorial  Report)  by  David  W. 
Vhitcembe,  Space  Technology  Laboratories,  12  April  1959. 
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Required J^ocity  -  ThlB  concept  Is  based  on  the  fact  that  at  each  space- 
tine  pdlnt  In  the  povered  flight  region,  a  required  velocity  vector 

Vp  -  ys,  t)  (29) 

aay  be  defined  and  computed  such  that  the  resulting  free  flight  trajectory  vill 
satisfy  four  general  guidance  constraints. 


Consider  for  exaaple  the  ICBH  trajectories  shovn  in  Figure  6.  Although  A1  and 
A2  arrive  at  point  A  at  different  times,  all  three  trajectories  can  be  made 
to  "hit''  the  target.  That  is,  after  a  period  of  time,  t^^  (the  time  of  free 
fli^t)f  the  three  trajectories  satisfy  the  three  guidance  constraints  that  R 
be  equal  to  .  Consider,  hovever,  the  two  trajectories  originating  at 


Earth 


Fig^e  6  THREE  FREE  FLICTT  TRAJECTORIES,  WHICH 
ALL  HIT  THE  DESIRED  TAHOET 

point  A.  In  order  to  specify  the  required  velocity  at  A  ,  an  additional 
guidance  constraint  must  be  iBqjosed.  This  should  be  obvious  since  it  pre¬ 
viously  has  been  shovn  that  only  two  degrees  of  freedom  are  required  to 
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control  impact  and  naturally  has  three  coii5>onentB.  For  example,  it  la 
possible  to  hold  the  total  time  of  flight  from  launch  constant.  This  has 
the  effect  of  making  an  eax^  fixed  target  stand  still,  vhich  scnetimes  has 
certain  adyantages  from  a  computational  standpoint.  For  this  set  of  con¬ 
straints,  the  missile  vlU  reach  the  target  vith  an  arbitrary  velocity  vhich 
depends  on  the  pazticular  burnout  conditions.  Soiae  of  the  other  possible 
constraints  are 

1)  Burnout  velocity  magnitude 

Z)  Burnout  velocity  elevation  angle  \ 

3)  Any  component  of  bximout  velocity 
h)  Burnout  energy 

5)  Burnout  angular  skooentum  magnitude 

6)  Velocity- ma^ltttde' at  the  target 

7)  Velocity  elevation  angle  at  the  target 
6)  Any  coeg>onent  of  velocity  ao  tbe  target 
9)  Time  of  free  flight. 

Some  Exaiiglee  -  It  Is  not  necessary  to  include  all  three  target  position 
coordinates  as  constraints.  For  one  satellite  launch  problem  it  vas  found 
desirable  to  use  the  following  constraint: 

1)  X  -  x^ 

3)  y  - 

U)  constant  total  time  of  fll^t. 

Here  «  Xij*  ^  «nd  y^  are  all  allowed  to  vary.  In  the  case  of  satellites 

which  have  a  number  of  thrust  periods  separated  by  coast  periods,  it  Is  usually 
desiirable  to  use  different  sets  of  constraints  during  the  various  guidance  periods. 

The  generality  of  the  required  velocity  concept  Is  Illustrated  by  the  fact  that 
the  simple  case  of  steering  a  vehicle  to  a  particular  velocity  at  burnout  is 
Included.  Here  the  constraints  are 

1)  X  -  Xj 

2)  y  -  y^ 

3)  i  -  ’.T 

t„..  0 
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In  this  case  I  the  target  position  vhlch  reduces  to  the  burnout  position  viU  vary 
depending  on  the  propulsion  system,  vlnds,  etc. 

The  above  examples  have  Illustrated  that  all  essential  guidance  information  is 
Included  in  the  specification  of  the  required  velocity  vector.  If  the  missile 
is  steered  in  both  pitch  and  yaw,  it  la  theoretically  capable  of  achieving 
exactly  the  requjxe^  .velocity  when  thrust  is  terminated,  and  hence  capable  of 
satisfying  the  four  guidance  constraints.  Since  satisfying  four  constraints, 
actually  means  causing  three  of  the  constraints  to  occur  simultaneously  with 
the  natural  occurance  of  the  fourth,  it  is  reasonable  that  this  degree  of  per¬ 
formance  should  be  attainable  with  the  three  degree  of  control  available  in 
most  rocket  vehicles,  l.e.,  pitch  steering,  yaw  steering,  and  thrust  tcrmluation. 
As  was  pointed  out  in  an  earlier  paragraph,  it  is  possible  to  satisfy  all  seven 
constraints  if  the  rocket  engine  is  throttleable. 

Another  Approach  -  Since  only  three  constraints  arc  actually  required  to 
cause  a  ballistic  missile  to  inqiact  at  a  target,  it  is  sometimes  convenient  to 
work  in  terms  of  a  two  component  (x  and  y)  required  velocity  vector.  In  this 
ewe  t  is  allowed  to  be  arbitrary  and  Vj^  becomes  a  function  of  i  as  well 
as  position  and  time.  Hence 

-  Vp(R,  t,  i)  (30) 

Since  t  is  arbitrary  pitch  steering  is  not  required  for  guidance  purposes, 
it  may  be  used  to  satisfy  antenna  look  angle  constraints  or  may  be  eliminated 
altogether  in  the  interest  of  simplicity. 

Some  of  the  presently  used  methods  for  computing  required  velocity  are  described 
in  Section  3- 

2. 3  The  Control  Function 

The  basic  function  of  the  control  portion  of  a  rocket  guidance  system  is  to 
steer  the  vehicle  in  a  way  that  will  cause  the  actual  velocity  to  become  equal 
to  the  required  velocity.  When  this  occurs  thrust  Is  terminated  and  the 

^^hc  w:>rk:  of  D.  MaePherson  nt  Space  T**chnology  Labor atorlen  haa  shown  that  it  is 
often  rosslble  to  satisfy  more  than  four  constraints,  by  steering  in  n  particular 
way.  His  work  is  not  covered  here. 
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guidance  systems  plays  no  further  roil  unless  there  are  additional  guided  phases. 
Since  the  steering  loops  are  feedback  systema,  it  la  convenient  to  define  an 
error  signal,  the  moat  coanonly  used  one  being  given  by  the  expression 

9  =  V  -  ^  (31) 

g  R 

where  V  is  equal  to  R  ,  and  where  V  is  called  the  velocity- to-be -gained 

g 

or  alternatively  the  velocity- to-go.  It  is  the  function  of  the  steering  loop 

1  * 
to  drive  V  to  zero;  the  steering  loop  is  actvially  a  velocity  control  aysteo 
^  '  S 

with  Vj^  as  the  connnanded  input. 

The  Steering  Loops  -  A  block  diagram  of  a  high  quality  inertial  guidance 

system  for  rocket  vehicles  is  shown  in  Figure  7.  As  can  be  seen,  both  the 

navigation  loop  and  the  guidance  equations  are  Integral  parts  of  the  steeri^ 

loop.  The  velocity  control  lav  (or  function)  must  be  chosen  so  that  some  time 

will  alvaya  exist  where  all  components  of  V  are  equal  to  zero.  Since  the 

S 

system  is  non-linear  and  time-varying,  the  development  of  the  proper  function 
is  not  a  simple  task.  While  Jb  alii  Stic  guidance  systems  are  usually  thought  of 
as  velocity  control  systema,  as  stated  above,  an  outer  loop  dxjes  exist.  The 
predicted  values  of  the  three  or  four  target  variables  of  interest  are  con¬ 
stantly  being  compared  with  their  specified  values.  This  is,  of  course,  done 
implloltly  by  the  guidance  equations.  The  result  of  the  comparison  is  the 
command  signal  for  the  velocity  loop.  In  practice  V  changes  much  slover 
than  V  ,  so  for  stability  studies,  the  former  can  be  treated  as  a  forcing 
function. 

As  was  pointed  out  in  Section  1.2,  the  relatively  fragile  nature  of  long 
range  rocket  vehicles,  usually  requires  that  the  steering  be  divided  into 
two  phases;  (1)  an  atmosjfieric  phase  and  (2)  a  vacuum  phase.  During  the 
first  phase,  which  corresponds  to  the  region  where  aerodynamic  effects  are 
appreciably  and  where  the  nominal  trajectory  is  a  zero  lift,  or  gravity,  turn, 
only  a  rather  "gentle”  type  of  steering  can  be  tolerated.  Any  violent  maneuver¬ 
ing  vin  result  in  excessive  structural  loads  which  may  cause  the  vehicles  des¬ 
truction. 

♦The  term  control  system  is  frequently  taken  to  mean  attitude  control  system. 

In  this  chapter ,  'control"  has  a  more  general  waning.  When  attitude  control 
is  meant,  it  will  be  so  specified. 
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Steering  Signali  -  The  net  result  of  the  "gentle"  steering  requirement  Is 

that  V  Is  visually  not  an  acceptable  steering  signal  duri:ig  phase  1,  because 
S 

of  the  violent  maneuvering  vhlch  It  can  cause.  The  tvo  approaches  which  have 
been  used  noat  successfully  to  date  are 

(1)  Progranalng  pitch  and  yav  attitude  as  functions  of  time. 

(2)  Steering  the  missile  on  some  moderately  well  behaved  function 
such  as  one  which  cosinanda^a  given  velocity  profile,  e.g. 

null  y  with  yav  cuctrol  and  coesnand  a  functional  relationship 
between  ic  and  z  with  pitch  control. 

Any  scheme  used  during  the  atmospheric  phase  besides  being  "gentle"  should  heep 
the  vehicle  as  close  as  possible  to  the  noninal  trajectory,  even  In  the  phase 
ot  perturbations,  so  that  the  required  region  of  guidance  equation  validity 
is  sdnimized. 

lAirlng  the  second  or  vacuum  steering  phase,  yaw  attitude  or  attitude  rate  is 
commanded  with  sene  function  of  ^  Pitch  attitude  or  attitude  rate  may  be 
commanded  as  In  the  atmospheric  phase  or  by  some  function  ot  ,  depending 
on  the  number  of  guidance  constraints  to  be  satisfied.  In  any  case,  thrust  Is 

terminated  when  ^  xreaches,  zero  or  nearly  so.  The  steering  problem,  Including 

- 

a  number  of  presently  used  approaches,  is  discussed  more  completely  In  Section 
Regardless  of  the  system  used,  excessive  maneuvering  should  be  avoided,  clnce 
It  Is  Inefficient  from  a  fuel  standpoint. 

This  section  has  described  a  general  type  of  rocket  vehicle  inertial  guidance 
system.  It  was  .ihown  that  the  scheme  has  three  main  functions:  (1)  Navigation, 
(2)  Computation,  and  (3)  Control.  Such  systems  can  theoretically  achieve  almost 
arbitrary  accuracy,  limited  only  by  the  capacity  of  the  airborne  computer.  In 
practice,  of  course,  component  errors  place  a  llsdt  on  the  achievable  performance. 
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3-0  GUIDANCE  EQUATIONS 

Given  the  tliaej_and  t.-he  position  and  velocity  of  a  vehicle,  a  guidance  systea 
aust  be  capable  of  generating  signals  'which  are  a  measure  of  the  vehicles  present 
ability  to  accoapllsh  Its  Intended  mission.  'Riese  signals  are  in  turn  used  to 
Bodily  the  flight  path  and  to  terminate  thrust  in  a  way  which  guarantees  that 
all  guidance  constraints  will  be  satisfied.  The  generation  of  these  error  or 
status  signals  I  which  was  referred  to  as  the  computation  function  in  Section  2, 

Is  accomplished  through  the  use  of  guidance  equations.  Two  of  the  more  common 
types  of  guidance  equations  are  discussed  in  this  Section;  both  are  based  on  the 
required-velocity  concept. 

In  general  the  closed-fora  computation  of  required  velocity  for  the  case  of  a 
rotating,  oblate  earth  %rith  atmosphere  is  not  possible.  Approximate  calculations 
of  more  than  adequate  accuracy  can  be  accoBrpllshed,  however,  in  a  number  of  ways, 
two  of  the  more  useful  ones  being 

1)  Explicit  Guidance  Equations 
Z)  Delta  Guidance  Equations 

From  a  practical  standpoint,  the  methods  differ  in  two  Important  respects: 

(1)  complexity  of  the  In-fllght  cootputations  and  (2)  amount  of  targeting  or 
pre -computation  required.  Basically  one  can  be  traded  for  the  other.  Explicit 
€»quations  are  rather  complicated  from  a  mechanisation  standpoint ,  but  require 
a  sdnifflua  of  pre-computation.  For  delta,  essentially  the  reverse  Is  true. 

The  selection  for  any  actual  system  ■will  depend  mainly  on  the  relative  importance 
of  these  two  points. 

3 •  1  Ex.plicit  Guidance  Equations 

The  mort  straight  forward  approach  to  the  guidance  of  a  biillstic  vehicle  Is 
to  work  directly  with  the  free  flight  equations  of  motion.  For  purposes  of 
Illustration  consider  again  the  case  of  a  laissile  traveling  in  a  vacuum  under 
the  influence  of  a  constant  gravity  field  (uiually  called  a  "flat  earth").  The 


STL/ni-60-0000-<2ll86 
PiM^e  25 


equations  of  ■otlon  as  prenrlously  given  in  Section  1  are 


X  «  X  +  X  t^^ 

o  O  ff 


'  -  ‘o  *  «  ‘ff 


(32) 

(33) 


2 

where  x  is  downrsnge,  t  is  up,  and  g  is  equal  to  32.2  fps  .  The  geometry 
Is  Bhovn  in  Figure  8.  Motions  out  of  the  trajectory  plane  are  not  considered 
for  the  present. 


Figure  8  -  FLAT  KARTH  EJCFLiaT  OUIBAIfCE  EDCAMPLE 


If  the  total  tine  of  flight  is  constrained  to  be  equal  to  T  then  the  x 

and  z  coaponents  of  the  required  velocity,  V_  ,  are  given  by 

^ —  n 


+  1/2  g(T  -  t^) 


vhere 


'ff 


-  T  -  t. 


W 

(35) 

(36) 


If  the  Blssile  is  steered  in  pitch  so  that  T  Is  equal  to  V-  and  if 

X  nX 

the  thrust  is  terminated  equal  to  V_  ,  then  the  aissile  vlll 

X  RX 

pass  through  point  (x^  ^  time  T 
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If  impact  time  is  of  no  interest  as  is  sometimes  the  case,  the  tir  of  free 
flight,  t^  ,  can  be  found  from  (33)  aa  a  function  of  x  and  j  as  given  by 


g 


(37) 


Substituting  the  expression  for  t 


ff 


for  V, 


Rx 


(37),  back  into  (32)  glvep  the  expression 


g(x^ 
■77? 


i  -  2g(t^  .  z^) 


(38) 


The  cojsputatlon  of  Vj^  Is  more  coeipllcated  than  In  the  constant  ‘ime  of  fll^t 
case,  but  the  need  for  pitch  steering  has  been  eliminated,  and  It  s  nov  only 
necessary  to  cut-off  the  engine  when  is  equal  to  V_  For  either  of  these 
simple  systems,  which  actually  approximate  the  case  of  a  short  rar -*e  missile, 
lateral  motions  can  be  handled  by  constraining  the  missile  to  fly  in  the  launch- 
target  plane,  l.e.,  by  nulling  y  and  y  by  means  of  yaw  steerir:".  The  case  of 
a  moving  (but  not  accelerating)  target  Is  easily  handled  by  the  ur  of  a  target- 
fixed  coordinate  system.  If  the  target  Is  accelerated,  but  in  a  Vnown  way,  the 
problem  can  still  be  solved,  but  the  equations  are  more  complex. 


"flie  Spherical  Earth  Case  -  The  more  important  and  considerab  more  com¬ 
plicated  case  of  a  missile  traveling  about  a  spherical  earth  will  row  be  Investi¬ 
gated.  The  effects  of  earth's  oblateness  and  the  re-entry  atmosph-re  (in  the 
case  of  a  baHlstic  missile)  are  relatively  small  effects,  which  rn  be  handled 
separately,  and  which  will  be  discussed  later.  If  the  trajectory  n  again  con¬ 
sidered  planar,  the  Lagranglan  function  for  the  missile  In  free  fi ight  Is  equal 
to 


L 


2^2  ^ 
r  0  ) 


+  0*^  ^ 


r 


(39) 
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vhere  r  and  0  are  earth  centered  Inertial  (ECI)  coordinates.  From  (39) 
the  equations  of  motion  are  found  to  be 


re 


GM 

r 


Equation  (Ul)  say  be  integrated  directly  to  give 


0 


■j 

r 


If  the  variable  p  «  r  is  introduced,  then 


dp  dr  ^ 


dt 


dt  dr 


lE 

dr 


Substituting  (U2)  and  (U3)  Into  (UO)  gives 


dr 


li  + 

3  7 


vhlch  can  be  integrated  so  that 

2  <? 

a  .  Cj 


Re-arranging  (4^)  and  taking  the  square  root  of  both  sides  leads  to 

T)  - 


2QM 

r 


i 

r 


But  since 


r  ■ 


dr  de 
de  dt 


dr  1 

de  “7 
r 


(40) 


(42) 


(U3) 


(U4) 


(45) 


(46) 


(U7) 
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equation  (1*6)  becomes 


(U8) 


whonc  integral  le 


(**9) 


Since  9  is  sirbitrary  it  can  be  set  equal  to  rero.  The  equation  can  then 
o 

be  re-written  so  that 


r 


(50) 


This,  however,  le  identical  in  form  to  the  equation  for  an  eUlpse  with  one 
focus  at  the  origin 

.  °(i  -  ch.  (51) 

^  1  -  e  slnO 


where  a  is  the  eemi-major  axis  and  e  Is  the  eccentricity.  By  equating 
similar  coefficients  in  (50)  and  (5l)»  constants  of  integration  are  related 
to  the  elliptical  constants  as  follows; 

C"  ^  YoTotd  -  J)  '52) 

c  .  .21  (53) 

2a 
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Therefore  using  (U2)  and  (53)»  cqxiation  (45)  can  be  reduced  to 

r^  +  {erf  -  .  GM(|  -  i)  (5**) 

» 

This  is  the  so  called  "vis  viva"  Integral  of  celestial  mechanics  vhlch 
■imply  shovs  that  the  speed  of  a  body  in  orbit  depends  only  on  its  distance 
from  the  center  of  force.  It  nov  becomes  a  problem  in  analytic  geometry  to 
derive  an  expression  for  the  velocity  required  to  hit  a  desired  inpact  point, 
since  (54)  clearly  shovs  it  is  necessary  to  have  only  the  semi -major  axis  of\ 
the  free  flight  ellipse  Joining  vehicle  and  target. 

The  Required  Velocity  Expression  -  In  order  to  derive  the  expression  for 
the  semi -major  axis  of  the-«lltpse,  it  is  necessary  to  use  the  following  two 
properties  of  sn  ellipse: 

(1)  The  s\im  of  the  distances  from  sny  point  on  the  ellipse  to  the  two 
foci  is  equal  to  the  major  axis  of  the  ellipse. 

(2)  A  line  normal  to  an  ellipse  bisects  the  angle  between  the  two  lines 
to  the  foci. 

The  geometry  of  the  situation  is  shown  in  Figure  9»  where  f  and  r^  are  the 
missile  snd  target  vectors,  s  is  the  range  angle,  and  y  is  the  angle  of 
with  respect  to  the  local  vertical.  From  the  figure  it  can  be  seen  that 


2  2 
x  +  y 

2 

=  ( 2ct  -  r^) - — 

(55) 

y 

*  r  +  (aa  -  r)coa2r  -  r  cos# 

(56) 

X 

«  (aa  -  r)  sln2r  - 

(57) 

•Slnilar  derivations  of  tbe  %is  viva"  Integral  can  be  found  In  almost  any 
book  on  classical  mechanics,  e.g. ,  McCuskey,  S.  V. ,  "An  Introduction  to  Advanced 
Dynamics",  Addlson-Wesley,  1959- 
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SubBtitutlon  of  (56)  into  (5^)>  th*  viva”  Integral,  gives 

.  221  - i_l_££5Lf -  (59) 

^  **  ”  eosSt)  -  cos*  +  coB(2ir  -  •) 

t 

vhich  is  an  expression  for  the  reqvilred  velocity  In  terms  of  the  missile  and 
target  distances  from  the  center  of  the  earth,  the  range  angle,  and  the  velocity 
vector  angle.  It  may  he  computationally  convenient  to  mechanize  (^)  as  a  vec¬ 
tor  equation,  replacing  cosines  vith  dot  products.  Thrust  termination  can  be 
comnanded  according  to  either  of  tvo  criteria 


(» 


>  (60) 


-V 


vfaere  and  V  are  both-scalar  quantities.  Equation  (60)  has  soam  advantages 
since  it  tends  to  be  linear  in  the  region  of  cutoff. 


Effect  of  Earth's  Rotation  -  Thus  for  an  equation,  (59)»  has  been  derived 
vhich  gives  the  velocity  required  to  hit  a  particular  point  ^  inertial  space. 
Most  targets,  hovever,  are  not  fixed  in  Inertial  space,  they  are  usually  either 
in  orbits  of  their  ovn,  or  they  are  fixed  to  the  earth  and  rotate  with  it.  In 
particular  the  ECI  rectangular  coordinates  of  an  earth  fixed  target  at  sosie 
future  time  are  given  by 


*t  *  *’t  ^t  *  ®,  ’■ff)  (62) 

CO.  9^  T„)  (63) 

•  r^  sin  (6U) 

vhere  0^  and  are  the  latitude  and  reference  longitude  of  the  target  and 
vhere  is  earth's  rate  and  T^j  is  the  time  of  free  flight. 
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The  above  equations  clearly  chov  that  in  order  to  predict  the  position  of  the 
target  at  impact,  the  time  of  free  flight  nust  be  known.  The  time  of  free 
flight,  however,  depends  on  V_  which  in  turn  is  a  function  of  the  Impact 

n 

position.  Hence,  an  Iterative  procedure  results  as  follows; 

1)  A  future  target  position  Is  assumed. 

2)  The  corresponding  required  velocity  is  computed. 

3)  The  elements  of  the  resulting  ellipse  are  computed. 

U)  The  time  of  flight  Is  computed. 

5)  A  new  futxire  target  position  is  computed. 

0)  The  pfocadure  is  repeated. 

Because  the  earth  moves  rather  slowly  (1000  fps  at  1+5°  latitude),  the  time  of 
flight  calculation  is  not  very  critical  and  the  Iteration  can  be  carried  out  on 
a  cycle  by  cycle  basis.  This  Is  not  necesseurily  the  case  for  targets  in  general. 
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It  Is  now  convenient  to  Introduce  another  frequently  used  orbital  paraneter  -  B  » 
the  eccentric  anomaly.  Tbe  geometry  which  defines  the  an^e  Is  shown  in  Figure  10. 


►  X 


Figure  10  -  ORBITAL  GEOMETRY  SHCWI50  MEAN  ANCKALY  ANGLE 

Concentric  circles  of  radius  o  and  b  are  drawn.  If  a  radius  vector  is 
drawn  at  angle  E  with  respect  to  the  x  axis,  and  if  parallels  to  the  x 
and  y  axis  are  drawn  through  the  points  where  the  radius  vector  intersects 
the  small  and  large  circles,  then  the  lnter>section  of  the  two  parallels  will 
be  a  point  on  an  ellipse.  The  eccentric  anomaly  is  a  way  of  designating  a 
particular  point  on  the  ellipse  Just  as  ia  9 

From  the  geometry  of  Figure^lOr  it  is  easily  shown  that 


r  - 


a  (1  -  e  cosE) 


(69) 
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or  that 

e  coaE  «  ^  "  3 

Tlierefore  (68)  becomes,  after  taking  a  square  root, 

_  ♦  ____ 

—  -  ^/aS  e  Bins  (71) 

Va 

Multiplying  (69)  by  its  time  derivative  gives 

rf  ■  e  8lnE(l  -  e  cosE)  E  (72) 

which  can  be  conbined  with  (71)  to  give 

•yS  -  E  (1  -  e  cosE)  (73) 

which  in  turn  can  be  integrated  directly  to  give  Keplers  equation 

-  e  sinE)  «  •JBi  t  C  (7U) 

where  C  is  the  constant  of  integration.  The  time  of  flight  to  ,the  target 
follow*  directly  from  (74)  and  is  given  by 

-  (GM)'  '  E)  -  e(BlnE^  -  slnE))  (75) 

where  E  and  E^  are  the  eccentric  anomolles  of  the  missile  and  target 
respectively.  Since  all  terms  of  (75)  are  easily  computed  from  orbital 
relationships,  the  inertial  position  of  the  target  can  be  obtained. 

Other  Effects  -  The  scheme  described  here  ia  designed  to  hit  a  target 
moving  in  a  known  manner,  such  as  a  point  on  the  earth  or  a  catellite.  It 
does  not  require  control  of  radial  velocity,  i.e.,  the  pitch  attitude  program 
can  be  arbitrary.  If  a  more  general  gtiidance  problem  such  as  the  placing  of 
a  satellite  in  orbit  is  to  be  solved,  then  radial  velocity  must  be  controlled 
and  the  equations  become  more  conplicated. 
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Thus  far  no  mention  has  been  made  of  ways  to  compensate  for  re-entry  and  oblate¬ 
ness  effects,  nor  have  vays  to  handle  Taotlons  out  of  the  orbital  plane  been  con¬ 
sidered.  The  first  two  effects,  It  turns  out,  are  relatively  small  and  can 
usually  be  handled  by  tabulated  target  offsets  which  are  functions  of  range, 

azimuth,  latitude,  etc.  Ihe  last  effect  is  more  complicated  in  that  some  steering 

\ 

is  necessary.  Beisically  rJJ.  that  is  required  Is  that  the  vehicle  have  no  velocity 
normal  to  the  plane  containing  the  vehicle,  the  target  at  the  time  of  impact,  and 
the  center  of  the  earth.  The  normal  velocity  is  given  by 

-  (f  X  r^)  .  V  (76) 


which  can  be  used  as  the  error  signal  for  yaw  ateeriiig  and  is  analogous  to 


Delta  Guidance  Equations 


Until  very  recently,  the  complexity  of  explicit  equations  made  them  very  unattrac¬ 
tive  for  use  in  inertial  guidance  systems.  Airborne  computers  of  rewonable  site 
simply  were  not  fast  enough  to  solve  the  equations  in  real  time,  e.g. ,  once  or 
twice  per  second.  It  was  therefore  necessary  to  find  other  ways  to  handle  the 
problem. 


Since,  as  has  already  been  stated,  vehicle  perturbations  are  relatively  small 
(in  the  neighborhood  of  a  few  percent)  it  was  only  reasonable  to  think  in  terms 
of  guidance  equations  which  are  power  series  expansions  about  a  nominal  trajectory. 

General  Development  -  A  power  or  Taylor  series  in  one  variable  is  of  the  form 

f(x)  *  f(a)  +  (x  -  a)  +  (x  -  a)^  ♦ -  (77) 

*  !•  the  variable  and  a  is  the  point  about  which  the  expansion  occurs. 

For  a  function  of  two  variables,  f(x  ,  y)  ,  the  series  is  of  the  form* 

- 

f(x,y)  -  f(a,b)  ♦  (x-a)  ^  +  (y-b)  ^ 

r  ? 

^  ^  ^  j  f(x»y)  ♦  —  (T8) 

a,b 

•See  any  advanced  calculus  book,  e.g.,  C.  R.  Wylie,  "Advanced.  Engineering  Mathematics," 
pp  596-bOU,  McGraw-Hill,  1951. 
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vtaere  x  asd  y  are  the  variables  and  a,b  Is  the  point  about  vhlch  the  series 
Is  expanded.  The  similarity  to  the  one  variable  case  Is  unmistakable.  The  ex¬ 
tension  to  three  or  siore  variables  Is  equally  clear. 


As  has  been  shovn  previous ,  V-  has  either  tvo  or  three  components  depending  on 

n 

the  number  of  guidance  constraints  to  be  satisfied.  In  the  more  general  three 
cosqponent  case,  each  component  Is  a  function  of  four  variables  -  x,  y,  z,  and  t. 
They  are,  of  course.  Implicit  functions  of  the  guidance  constraints  themselves  also. 


If  expansions  are  to  be  found,  then  an  expansion  point  must  be  selected.  The  first 
approach  sdght  be  to  select  points  all  along  tbe  nominal  traO^ctory  and  to  program 
them  as  a  function  of  tlste.  Coefficients  vould  have  to  be  determined  correspond¬ 
ing  to  each  expansion  point  and  the  three  expansions  (Vj^  ,  and  vould 

be  time  varying.  If  the  number  of  expansion  points  used  was  very  large,  the  air¬ 
borne  computer  storage  requirements  could  easily  become  excessive. 


A  little  reflection  on  the  problem,  however,  soon  leads  one  to  tbe  conclusion 
that  the  expansions  need  to  be  highly  accurate  only  In  the  Immediate  vicinity 
of  burnout.  Hence,  only  a  single  point  and  a  single  set  of  expansions  need  be 
used,  at  least  for  any  one  guidance  phase.  Since  the  nominal  burnout  point 
(x^  >  y^  >  Sq  »  t^)  la  usually  coasldered  to  be  the  most  likely  burnout  point, 
(or  vector)  It  Is  usually  selected  as  the  expansion  p)olnt.  Three  expansions 
alsdlar  to  the  following  Vj^  expression  therefore  result: 


'Rx 


*0  ^ 


k^i^  +  k  dy 
XX  xy^ 


k  ^ 
xz 


‘'xt^^ 


k  £ix 

XXX 


k  £ix£iy  *  k  ♦ .  (79) 

xxy  ^  xyy  ' 

where  ^  ■  (x-x^),  ^  ■  (y-y^),  etc.,  the  delta  quantities,  of  course,  giving 
the  equations  their  name.  The  coefficients,  k^  ,  k^  .  etc.,  are  basically 
partial  derivatives  and  are  defined  by  tbe  corresponding  terms  In  (76).  Usually 
linear  and  some  second  order  terms  are  required  depending  on  the  probable  size 
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of  the  "burnout  box"  and  the  accuracy  required.  The  number  of  temuj  required 

I 

will  depend  also  on  the  alas  coefficients  in  the  direction  of  the  expansions 
validity. 

A  Flat  Barth  Exaaple  -  In  order  to  Illustrate  the  working  of  delta  equa¬ 
tions,  a  flat-earth,  two-dimensional  case  will  again  be  considered.  The  nominal 
trajectory  parameters  are  those  of  the  V-2  missile  from  Section  1  end  are  as 
listed  in  Table  3*  , 


Table  3  -  V-2  TWUKCTOPT  PARAMETERS 


Parameter 

Value 

Horizontal  BO  Position 

85,000  feet 

Vertical  BO  Position 

*o 

123,700  feet 

Horizontal  BO  Velocity 

« 

3,500  fps 

Vertical  BO  Velocity 

*o 

4,800  fps 

Time  at  BO 

t 

o 

70  seconds 

Horizontal  Target  Position 

1,212,000  feet 

Vertical  Target  Position 

“t 

0  feet 

Time  at  Impact 

T 

392  seconds 

The  constant  time  of  flight  constraint  will  be  employed  (T  -  392)  and  both 
linear  and  quadratic  tennB  will  be  used.  The  coefficients  can  be  obtained  by 
performing  the  appropriate  partial  differentiations  on  the  explicit  expressions 
for  Vj^  and  Vp^  ,  (3*^)  and  (35)*  The  coefficients  are  suaiBarlxed  In 
Table  4.  It  can  be  seen  that  while  there  are  l8  possible  linear  and  quadratic 
termr  in  the  two  expansions,  only  8  are  non-xero  for  the  flat-earth,  constant¬ 
time  -of  -flight  case.  The  biggest  reduction  la  due  to  the  absence  of  cross 
coupling  between  the  x  and  x  channels. 
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Table  U  -  LIWEAR  AND  QUADRATIC  DELTA  COEFFICIENTS  FOTl  V-2  TRAJECTORY 

(Constant  Tiae  of  FliRht) 
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Usln^  the  V-2  paraffleters^  the  V_  and  V  expansions  become 

Rx  Ri 

^Rx  “  ‘  -  85,000)  +  10.9  (t  -  70)  (80) 

-  9.6  X  10*^(x  -  85,000)(t  -  70)  4.  3.U  X  10*'^  (t  -  70)^ 

*  U800  -  3.1  X  10"^  (t  -  123,700)  .  17.3  (t  -  70)  (81) 

-6  p 

-  9.6  X  10  (z  -  123,700)(t  -  70)  -  3-7  X  10"  (t  -  70) 


Vlthln  the  airborne  computer,  It  Is  usually  convenient  to  group  the  terma  some¬ 
what  differently,  lunging  the  constants  all  into  a  single  tern. 

In  order  to  give  aone  idea  of  the  accuracy  of  the  equations  as  well  aa  the 
relative  lBq>ortance  of  the  vtu'lous  tema,  the  equations  have  been  evaluated 
for  three  tats  of  perturbations,  one  of  which  corresponds  to  the  launch  point. 
Thr  results  are  suflBMrized  in  Table  5.  The  launch  point,  of  covirse,  is  not 
an  expected  burnout  point  and  is  only  included  for  the  purpose  of  placing 
some  sort  Of  an  upper  bound  on  equation  error.  For  the  other  randomly  selected 
perturbations  (which  are  abnormally  large  by  any  standards),  the  equations  per¬ 
form  quite  well.  Even  with  linear  terms  alone  the  V-2  impact  error  would  be 
well  under  a  half  mile  if  the  missile  were  to  bum  out  at  these  points.  The 
example  also  graphically  Illustrates  the  relative  importance  of  linear  and 
higher  order  terma. 

The  More  C>eneral  Caae*"^  While  the  flat -earth  example  Illustrates  very 
well  the  functioning  of  delta  equations,  it  gives  a  somewhat  over-simplified 
picture  of  the  effort  that  must  go  into  obtaining  the  coefficients  -  a  Job 
usually  known  as  "targeting".  For  the  flat-earth  case,  simple  formulas  are 
available  for  the  required  velocity  and  the  partial  derivatives  are  rather 
quickly  obtained,  nie  Vp  formulas  themselves  are  simple  enough  to  pro¬ 
bably  obviate  delta  equations. 

♦To  tHe  writer's  knowledge,  the  first  successful  application  of  delta  guidance 
equntii,'ns  to  long  range  rockets  was  a  result  of  the  Joint  efforts  of  J.  Caroll 
of  American  Bosch  Arma  Corporation  and  F.  Baskin  and  T.  W.  Layton  of  Space 
Techno]  Ofty  Laboratories. 
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Table  5 

BEHAVIOR  OF  DELTA  GimiANCE  EQUATIONS  FOR  V-2  TRAJECTORI 
(ALL  VELOCITIES  IN  FPS;  CALCULATIONS  TO  NFJU^EST  1  FPS) 


I  Launch 

Ax  -  20,000  ft  1  Ax  -  -  10,000  ft  Ax  -  -  85,000  ft  | 

Ax  -  30,000  ft  i  Ax  •  -  15,000  ft  Ax  «  -  123,700  ft 


At  -  - 

10  sec 

■  At  ■  +  10  sec 

At  «  -  70  sec 

^Rx 

^Rx 

^Rx 

^Rx 

_ _ _ ^ 

3500 

4800 

3500 

4800 

3500 

-  62 

0 

31 

0 

263 

0 

0 

-  93 

0 

u6 

Ci 

383 

-109 

173 

109 

-173 

-763 

12U 

2 

0 

1 

0 

-  57 

0 

0 

3 

0 

1 

0 

-  83 

3 

0 

3 

0 

l66 

-  18 

3329 

4880 

3640 

4674 

j  3000 

3334 

4883 

3644 

4674 

3109 

6293 

3333 

4882 

3644 

4675 

3092 

6311 

-  4 

-  2 

-  4 

-  1 

-  92 

83  ' 

.  J 

>  1 

*  1 

0 

0 

*  17 

-  18  i 

1 

_ i 
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When  a  rotating,  oblate  earth  vith  re-entry  ataoaphere  la  considered,  the  picture 
changes.  While  explicit  equations  for  a  spherical  earth  vithout  atmosphere  are 
presented  in  Section  3-l>  it  should  be  pointed  out  that  even  these  approximate 
equations  are  iterative  in  nature  and  can  not  be  readily  differentiated  in  order 
to  obtain  delta  coefficients.  The  procedure  which  has  evolved  to  generate  delta 
constants  basically  consists  of  the  following  two  steps; 

(1)  The  required  velocity  (V„  ,  V  ,  and  V_  )  corresponding  to  each  of  a 

KX  Hy  RZ 

large  number  of  x  ,  y  ,  o  ,  and  t  perturbations  is  found  by  using 
iterative  procedures  on  a  digital  coaiputer  free  flight  simulation.  The 
simulation  is  a  veiry  precise  model  of  a  rotating,  oblate  earth  with  at¬ 
mosphere.  The  points  are  selected  in  a  way  which  makes  them  correspond 
as  closely  as  possible  to  actual,  realltable  powered  flight  burnout 
points.  The  will,  of  course,  be  a  function  of  the  constraints 

used. 

(2)  The  empirically  generated  V_  data  and  the  corresponding  perturbations 

K 

are  used  to  generate  the  delta  coefficients  by  means  of  least  square 
fitting  procedures . 

The  procedure  for  generating  each  of  the  three  expansions  Is  the  same  as  the  one 
for  V  which  is  summarized -by  " 

KX 


_(h)N 
PU  ) 


AV, 


Rx 


(k) 


2 

-♦Min  (82) 


where 


Rx 


number  of  data  points 

•  •  •  •  and  At  from  the  kth  data  point, 

functions  of  Pj^)  like  Ax,  Ay,  AxCy,  etc. 

coefficients  of  AJ  In  the  expansion 

number  ofTlelta  terms  in  AlTp 


The  functions,  AJ  ,  and  the  number  of  terms,  n  ,  are  chosen  to  obtain  the 
required  degree  of  accuracy.  It  usually  turns  out  that  n  is  Boine%rtiat  less  for 
'^Ry  '"^Rz  fof  Vpjj  •  The  actual  differences  depend  mainly  on 

the  coordinate  system  and  the  constraints  employed. 
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The  linear  V  coefficients  can  also  be  obtained  from  the  linear  target  miss 
R 

expansions.  Since  these  expansions  are  more  easily  obtained  (no  free  flight 
iterations  are  required),  the  procedure  is  sometimes  used  vhen  only  the  linear 
coefficients  are  required  or  vhen  it  is  desirable  to  simplify  the  fitting  pro¬ 
cess.  Consider  for  simplicity  the  tvo  dimensional  case.  The  linear  miss 
expansions  are  of  the  form 

dM  SM 


an  SM  hH 

vhere  M  and  M  are  the  x  and  t  deviations  from  x  ,  z.  at  time 

X  S  ^  V 

equal  to  T  (or  as  any  other  third  constraint  Is  satisfied).  If  the  misses, 

fi  ,  are  set  to  zero,  then  the  resulting  tvo  linear  equations  can  be  solved 

slBBiltaneously  for  Cx  and  /^z  The  resulting  expansions  in  ^  ,  and 

^  vlll  be  Identically  the  same  as  the  linear  V-  and  V_  expansions.  As 

RX  R£ 

an  example,  the  £ix  coefficient  In  the  V-  expansion  Is 

RX 


\;5rw 
(  sr  sr 


dM  SM  X 

x  z  \ 

5ir^; 

dM  dM 

z  X 

W  51”.' 


(85) 


^e  results  are  similar  in  the  three-dimensional  case,  except  that  numerator 
and  denominator  contain  six  triple  products,  rather  than  two  double  products. 
Numerically,  the  computations  are  trivial  in  either  case. 


Summary  -  Explicit  and  delta  guidance  equations  have  been  described  and 
examples  of  each  have  been  given.  While  these  are  not  the  only  guidance  equa¬ 
tions  currently  being  used,  they  do  typify  the  two  larger  classes  into  which 
all  guidance  equations  might  be  grouped,  namely  total  equations  and  perturba¬ 
tion  equations.  Total  equations  require  less  pre-computation  and  will  work 
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for  larger  vehicle  .and  environmental  perturbations.  Perturbation  equations, 
hovWer,  are  simpler  to  mechanize  and  are  Inherently  more  flexible,  sinca  all 
guidance  constrainta  can  be  changed  simply  by  changing  the  constants.  As  has 
been  stated  previously  the  choice  vlll  depend  mainly  on  operational  requirements. 
ITrequently  It  reduces  aLnost  to  a  matter  of  personal  choice. 
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U.O  VEHICLE  STEERING 

In  the  previous  section,  It  vas  shown  that  nethods  exist  for  continuously 

computing  the  required  velocity  for  a  rocket  In  flight.  The  third  and  final 

function  of  the  guidance  system,  as  described  In  Section  2,  is  the  control 

function,  i.e.,  the  vehicle's  flight  path  must  be  directed  in  a  way  which  will 

cause  all  three  components  of  the  veloclty-to-be-gained  (V  *  V_  -  V)  to  reach 

g  n 

zero  slnultaneously..  For  reasons  developed  in  Section  2.3>  the  steering  is 
usually  broken  dovn  into  atmospheric  and  vacuum  phases.  In  this  section,  the 
steering  problem  is  discussed  In  some  detail  and  a  few  of  the  currently  used 
methods  for  each  phase  are  described. 

U.l  The  Atmospheric  Phase 

During  the  atmospheric  pheise,  the  steering  of  a  rocket  vehicle  is  basically 

"open  loop",  I.e.,  V  is  not  used  expllclty  to  control  the  flight  path. 

S 

Instead  some  less  violent  variable  such  as  attitude  or  velocity  is  used  to 
maneuver  the  vehicle  in  a  way  which  will  keep  it  as  close  as  possible  to  the 
nominal  trajectory  without  causing  excessive  structural  loading. 

The  starting  point  for  most  atmospheric  steering  schemes  is  the  open  loop  or 
reference  trajectory,  usually  a  "kick"  trajectory.  This  non-physical  trajectory 
is  generated  on  a  digital  computer  by  causing  the  missile  to  rise  vertically  for 
some  period  of  time,  at  the  end  of  which  the  vehicle's  attitude  and  velocity 
vector  are  instantaneously  rotated  downward  by  an  amovmt  known  as  the  "kick" 
angle.  From  this  point  on  until  the  missile  is  essentially  out  of  the  atmosphere, 
the  thrust  vector  is  caused  to  be  directed  along  the  velocity  vector  and  the 
missile  flies  a  gravity  turn.  Beyond  the  atmosphere  some  arbitrary  attitude 
progrma,  usually  a  constant  angle,  is  flown. 

Since  the  velocity  vector  of  a  physical  rocket  u.-nnot  bo  Instantaneously 
"kicked",  the  transition  from  vertical  rise  to  gravity  turn  actually  takes  place 
over  a  period  of  time  known  an  the  transition  turn,  the  length  of  which  la 


STL/TJf-60-0000-<3Wl86 
P««e  *^5 


aeaBured  in  tens  of  seconds.  An  anpie  of  attack  mat  of  necessity  exist  during 
this  period,  but  fortunately  this  occurs  before  the  dynamic  pressure  has  not 
had  time  to  build  up. 

The  Bimplest  approach  to  steering  in  the  atmosphere  is  to  program  vehicle 
attitude  or  attitude  rate  as  a  function  of  time.  Yav  attitude  is  held  to  xero, 
and  pitch  attitude  is  coonanded  as  it  occxirs  on  the  "kick"  trajectory.  The 
"kick"  angle  is  actually  caused  to  occur  slowly  over  a  period  of  a  few  seconds. 
A  typical  pitch  rate  program  is  shown  in  Figure  11.  The  program,  which  In  this 
case  changes  in  steps  can  be  generated  mechanically  or  electronically.  The 


e 


c 

(deg/sec) 


Time  in  Seconds 


Figure  11  -  TYPICAL  PITCH  RATH  PROGRAM 

corresponding  pitch  attitude  program  could  also  be  used.  In  that  case  a  more 
continuous  program  is  required,  since  attitude  is  conaaanded  directly,  without 
the  smoothing  effect  of  an  integration. 
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Velocity  Steering  -  Soaievhat  better  performance  Is  obtained  by  steering 
the  mltalle  on  velocity.  The  velocity  may  or  may  not  Include  Integrated  gravita¬ 
tional  acceleration.  Consider  first  the  case  of  thrust  velocity  alone,  l.e., 
the  Integrated  output  of  acceleronetera.  The  coordinate  system  la  shown  In 
Wgure  12  where  x  y  t  (down  range,  left,  and  up)  Is  a  launch  fixed  Inertial 
coordinate  system.  If  the  x  axis  la  rotated  counter-clockwise  by  an  angle  \  , 


Figure  12  -  COORDINATE  STSTTM  FOR  THRUST  VELOCITY  STEERING 

the  u  axis  Is  obtained.  The  Integrated  outmit  of  an  accelerometer  mounted 
along  u  will  differ  from  u  at  any  Instant  of  time  by  the  Initial  value 
of  u  and  the  Integrated  gravity  term.  If  Is  equal  to  90°  minus  the 
constant  attitude  angle,  It  can  be  shown  that  for  typical  open  loop  trajectories 
the  actual  Integrated  accelerometer  output*,  u^  ,  will  have  a  tine  history 
of  the  fora  shown  In  Figure  13. 

•Aciualiy,  o^  course,  h  may  be  the  outputs  of  two  or  more  acceleronetera 
suitably  coetblned. 
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Figure  13  -  IirreORATED  ACCELERCKETTER  OOTPOT  FOR  TTPICAL  TRAJECTOTT 


The  shape  of  the  curve  suggests  an  exponential  approxlajitlon  of  the  form 

t 

m 

«  Uj.  (1  -  e  ^p)  (86) 

If  the  vehicle  can  he  caused  to  fly  such  that  (86)  is  realixed,  then  a  gravity 
turn  vlll  result.  This  can  he  acccopllahed  hy  comandlng  a  pitch  rate 


0 

c 


-  *  % 
P  p  A  A 


(87) 


If  l8  small,  then  (86)  is  an  exact  solution  of  (87) •  If  it  is  desirable 
to  command  attitude  rather  than  rate  then  (87)  can  be  integrated  and  the  command 
becomes 


e 

c 


K„[t  u  ♦  u  - 

p  p  a  a 


Ujt  + 


C] 


(86) 
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It  may  be  necessary  to  fit  the  ii  (t)  vs  t  curve  with  two  or  more  exponentials 
and  to  change  constants  at  some  convenient  point  such  sis  a  staging  point.  This 
is  especially  true  if  this  type  of  steering  is  to  be  used  over  the  entire  flight 
rather  than  for  the  vacuum  phase  o.  ly.  Velocity  steering  in  yaw  is  accomplished 
by  an  attitude  command  of  the  form 

♦e  -  -  'V'V*  * 

Here  there  is  no  forcing  function  and  the  yaw  steering  loop  is  actually  a  nulling 
loop. 

A  somewhat  laore  general  pitch  velocity  steering  scheme  consists  of  controlling 
the  flight  path  of  the  vehicle  so  that 

V  -  f.(V^  ,  X  ,  z  ,  t)  (90) 

Z  XX 

Again  a  vaxdety  of  coordinate  systems  with  aund  without  gravity  may  be  employed. 

An  esurth  fixed  coordinate  system  works  particularly  well,  since  the  earth's 
atmosphere  rotates  with  the  earth  and  therefore  is  at  rest  with  respect  to  the 
coordinate  system.  The  steering  necesseuy  to  realize  (90)  can  be  accomplished  by 
attitude  perturbations  given  by 

50^  -  f^(V^  ,  X  ,  t  ,  t)]  (91) 

where  60 /Kp  is  a  small  nximber.  At  the  same  time  the  nominal  pitch  program  can 
be  written  as  a  function  of  one  or  more  variables,  e.g. 

0  ■  f-(V  ,  X  ,  z  ,  t)  (92) 

nc  2'  X  ' 

Combining  (91)  and  (92)  to  get  the  total  pitch  attitude  angle  command  gives 


0  -  0  +  60  -  f  (V  ,  V  ,  X  ,  z  ,  t) 

c  nc  c  j  X  z 
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In  practice  it  la  frequently  found  that  the  x  ,  t  ,  and  t  terns  are  not  very 
effective  and  an  expansion  In  and  may  be  all  that  is  neee8saiTr»  «•£• 


0 


K 


X  V 

XX  X 


(91*) 


The  coefficients  are  obtained  by  selecting  a  nunber  of  points  from  the  standard 
trajectory  and  using  least  square  fitting  procedures  as  vas  done  for  delta  guidance 
constants.  The  difference  here  Is  that  It  Is  not  necessary  to  use  points  from 
perturbed  trajectories. 


Velocity  steering  as  typified  by  (88)  and  (9^),  although  more  complicated,  has 
at  least  tvo  advantages  vhen  compared  vlth  attltude>tlme  programming.  Plrst  the 
mlcsile  angle  of  attack  due  to  a  vind  snear  Is  reduced  sooevhat.  This  Is  due  to 
the  fact  that  the  dynamic  action  of  the  loop  causes  the  vehicle  to  "weathercock" 
into  the  wind.  The  second  and  primary  advantage  Is  that  the  trajectory  perturba¬ 
tions  due  to  vehicle  perturbations  are  greatly  reduced.  Vhen  any  type  of  pertur¬ 
bation  guidance  equations,  e.g.  delta,  are  used,  this  fact  results  In  a  reduction 
in  the  complexity  of  the  equations  and  hence  a  slaqiler  computer  mechanization. 

On  the  debit  side  is  the  fact  that  this  type  of  steering  tends  to  "work"  the 
attitude  control  system  harder  and  as  a  res\ilt  stability  margins  at  vibration 
frequencies  may  be  slightly  reduced. 


h.2  The  Vacuum  Phase 


Once  the  vehicle  is  out  of  the  atmosphere,  structural  constraints  can  be  relaxed 
and  the  steering  system  can  begin  to  perform  Its  primary  function,  fiat  of  re¬ 
ducing  V  to  zero.  A  steering  or  coaqputatlonal  coordinate  system  Is  usually 
cs 

selected  which  aligns  the  x  axis  more  or  less  vltL  the  desired  thrust  direction 

during  the  latter  portion  of  powered  flight.  The  signal  V  then  represents 

the  prlncii>al  component  of  veloclty-to -be -gained  and  V  and  V  are  con- 

gy  gz 

sidersbly  smaller. 


•The  coordinate  system  used  In  Section  3  had  the  s  axis  aligned  with  launch 
vertical  for  purposes  of  simplicity.  If  pitch  steering  Is  done.  It  Is  usually 
necessary  to  tip  the  coordinate  system. 


stl/tn-6o-oooo-<;ri86 
Page  50 


One  method  vhlch  will  cause  the  y  and  i  components  of  V  to  vanish  at 

6 

cutoff  Is  to  drive  these  quantities  continuously  tovard  tero  from  the  "initiation 

of  guidance"*.  This  can  be  accomplished  by  interpreting  V  and  V  as  error 

sy  K* 

signals  which  are  to  be  nulled  by  suitable  attitude  or  attitude  rate  commands  to 
the  autopilot  or  attitude  control  system.  Stability  considerations  usuadly  re¬ 
quire  that  some  type  or  rate  Information  also  be  Included  for  damping  p\irposes. 

A  typical  pitch  attitude  command  would  be  of  the  form 


e 


c 


K,V 
1  g* 


(95) 


where  9^  la  the  nominal  pitch  angle.  If  a  rate  attitude  command  is  to  be  used 
It  would  be  of  the  fom 


a)_ 


% 


(96) 


The  yaw  commands  would  be  the  same  only,  of  course,  V  would  be  replaced  by 

8* 

.  In  order  to  reduce  certain  steady  state  errors  It  is  often  desirable  to 
Include  also  an  integral  term  in  the  steering  expressions. 


An  Important  Modification  -  vniile  a  scheme  as  outlined  above  is  capable  of 
driving  f  to  xero  at  cutoff,  it  has  one  basic  shortcoming.  EVen  though  the 
y  and  x  components  of  V  are  small  compared  to  V  at  guidance  Initiation, 

o  g3C 

they  are  still  sufficiently  large  to  cause  excessive  pitch  and  yaw  coamiands  and 
hence  «ubstantlsl.'«aaeuverlng  of  the  vehicle.  As  a  result  even  a  nominal  vehicle 
Is  caused  to  deviate  considerably  from  the  reference  trajectory  and  this  is 
objectionable  for  two  reasons:  (1)  It  results  in  a  non-optimum  use  of  propellants 
end  (2)  It  means  that  the  required  velocity  computation  must  be  accurate  over  a 
larger  region. 


•'‘initiation  of  guidance"  is  a  term  often  used  to  Indicate  the  time  at  which 
the  vehicle  steering  Is  switched  from  the  atmospheric  phase  to  the  vacuum  phase. 
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The  reason  for  the  occurance  of  the  above  maneuvering  becomes  clear  vhen  It  Is 
realized  that  under  reference  conditions  the  y  and  x  components  of  V  differ 
substantially  from  xero  prior  to  cutoff.  Figure  lU  gives  the  general  shape  of  the 
three  coBiponents  for  a  typical  long  range  rocket  vehicle.  The  actual  shapes  will 


Velocity 


Time 

Figure  lU  -  OEWERAL  SHAPE  (F  V  COMPCWnTTS 

s 

depend  on  the  trajectory,  coordinate  system,  and  guidance  constraints.  The  net 
result  of  this  finding  is  that  not  V_  ard  V  ,  but  rathe’*  their  deviations 
fz*om  nominal  should  be  used  as  steering  error  signals.  If  this  is  done  non- 
noBiinal  missiles  still  "vander"  to  some  extent,  but  not  excessively  so. 


If  the  standard  values  of  V  and  V  are  to  be  progranaed,  the  question 
arises  as  to  the  Independent  variable  against  which  they  should  be  prograasied. 
Time  or  velocity  might  be  used,  but  further  consideration  leads  to  V  as  the 
best  choice.  The  pitch  and  yav  error  signals  are  then  of  the  form 


gx  gx 


f(V  ) 
'  gx' 


(97) 


(98) 
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where  the  ronction  aay  be  aiirply  V  nultlpilcd  by  n.  constant,  but  usually 
Includes  higher  order  tenaa.  The  function  will  be  equal  to  sero  at  cutoff; 
hence  the  erior  slgnalc  become  erual  to  V  and  V  at  cutoff.  This,  of 
course,  is  the  dominant  reaaon  for  using  as  the  programming  variable. 

CrosB  Product  Steering  -  The  steering  method  described  in  the  preceedit 


Cross  Product  Steering  -  The  steering  method  described  in  the  preceeding 
paragraphs  depends  on  the  exlstance  of  a  standard  or  reference  trajectory,  fre¬ 
quently,  especially  when  explicit  guidance  equations  are  used,  the  reference  tra¬ 
jectory  concept  loses  Its  meaning.  It  Is  then  that  more  general  steering  siethods 
must  be  employed.  The  best  known  of  these  is  cross  product  steering  where  the 

vehicle  is  caused  to  have  an  attitude  rate  proportional  to  V  x  V 

8  8 

The  fundaiiental  Justification  for  the  method  begins  with  the  definition  of  V 


Vg  -  -  V 


which  when  differentiated  with  respect  to  time  gives 


■  ''r  - 


(100) 


Since  V  is  adsslle  acceleration  which  during  powered  flight  is  due  to  thrust 
and  gravitational  accelerations,  (100)  can  be  rewritten 


Vg  .  ■-  0 


(101) 


where  I  is  the  unit  vector  in  the  direction  of  a^  ,  the  thrust  acceleration. 


Let  u  be  a  unit  vector  along  V  .  Therefore 

g 


(102) 


V  u  ♦  V  u 
8  g 


(103) 


•This  Section  is  based  mainly  on  work  by  J.  M.  Bachar,  F.  Baskin,  and  D.  W. 
Vhltcoabe  of  Space  Technology  Laboratories. 
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Since  u  la  a  unit  vector,  It  la  true  that 


(0 


u  X  u 


(110) 


and  hence 


U) 


V  V  V 

JL  X  -S-fi 

V  „  2 

*  ''g 


V  V  V 
*  "g 


(m) 


The  second  tern  la  Identically  lero  so  the  turning  rate  of  the  vector  is 
equal  to 


OQ 


V  X  V 
V 

g 


(112) 


If  the  thrust  vector  is  to  be  kept  aligned  with  the  V  vector  (C  •  u)  , 

D 

then  the  alsslle  ttimlng  rate  nust  be  equal  to  (112). 

An  examination  of  (112)  reveals  that  the  V  vector's  turning  rate  becomes 
infinite  at  the  instant  of  cutoff.  Since  no  physical  missile  hsis  this  capability, 
there  will  always  be  some  steering  error  at  cutoff  with  this  approach.  This  situa> 
tion  appears  to  be  analogous  to  a  pure  pursuit  course  where  infinite  turning  rates 
are  also  encountered  at  Ispact.  In  that  case  the  situation  is  alllviated  by  SKxlify* 
the  course  to  include  a  le$ul  angle.  In  the  case  of  steering  a  similar  effect 
is  achieved  by  bi_8ialng_out  the  nominal  error,  since  the  error  turns  out  to  be 
relatively  Independent  of  trajectory  perturbations. 

In  practice  the  cooraanded  vehicle  turning  rate  Is  of  the  form 


(o^  -  K[V^  X  V^) 

c  g  g 


(113) 


where  K  may  or  may  not  be  a  function  of  V  The  pitch  and  yaw  coBoanded 

g 

rates  will  be  equal  to 


®n  "  - 

P  _ _ _  -gX  gK  gt  gX 


u)  -  K(V  -  V  V  ] 

y  gx  gy  gy  gx 


(Ilk) 


(115) 
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The  corresponding  roll  rate  contains  redundant  Information  and  Is  not  used.  jUi 
equivalent  system  vhlch  coomands  angles  rather  than  rates  is  isosslble,  but  vlll 
not  be  descr^-bed. 

The  cross  product  steering  schem^causes  the  vehicle  to  fly  a  nearly  constant 
attitude  trajectory  'dilch  Is  desirable  from  a  fuel  standpoint.  At  guidance 
Initiations >  the  system  ccomanda  rather  large  rates  vtntll  the  proper  attitude 
Is  obtained  and  makes  only  small  changes  thereafter.  The  system  does  not  attempt 
to  remove  y  or  z  eoatponents  early >  rather  It  attempts  to  drive  all  components 
to  zero  simultaneously. 


ii.3  Thrust  Termination 


It  has  been  shown  that  It  Is  possible  to  steer  a  rocket  veh5.cle  In  a  way  which 

guarantees  that  all  three  components  of  the  velocity-to -be -gained  will  reach 

zero  at  the  same  instant  time  thrust  la  terminated.  If  the  steering 

system  Is  perfect,  l.e. ,  all  three  components  of  V  become  uniquely  zero, 

^  (S 

then  the  magnitude  of  V  or  the  magnitude  of  any  coaq>onent  of  V  Is  useful 

®  8 
as  s  cutoff  signal.  In  a  practical  system  where  steering  errors  will  occur  due 

to  both  static  and  dynamic  effects,  the  signal  la  usually  given  when 


8X 


<  K. 


CO 


(116) 


is  some  allowabla  tolerance  on  cutoff  velocity. 


where 

ponent  Is  chosen  for  a  nxiaber  of  reasons: 


The  V  com- 


(1)  The  coordinate  system  Is  usually  chosen  so  that  (V  |  Is  nearly  equal 

II  8*' 

(2) 

(3)  If  either  V _  or 


to  |V^|  . 

▼gjj  usually  already  available  since  it  Is  necessary  for  steering 
purposes. 


gy  -  were  used,  smaU  errors  In  these  quantities 

could  cause  large  cutoff  errors,  since  xisually  V  »  V  or  V 

8X  gy  gz 

further  consideration  of  the  cutoff  problem  depends  on  the  type  of  engine  et^iloyed. 
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3y«teii  VI th  Liquid  Rn^ln»»B  -  Timist  i*  temlnated  In  a  liquid  propellwJt 
crecinf*  by  closing  valves  in  the  fuel  and  oxlditer  lines.  Since  the  lines  pa^t 
the  valves  and  the  thrust  chairber  itself  contain  a  certain  amount  of  the  liquids, 
thrust  actually  continues  for  a  fraction  of  a  second  after  valve  closure,  a 
typical  value  being  0.2  seconds.  If  the  acceleration  near  burnout  is  10  g's,  and 
a  linear  decay  is  assumed,  32  fps  is  added  after  valve  closure  corresponding  to 
32  miles  miss  in  the  case  of  a  5500  mile  mlseile. 

A  known  residual  Impulse  can  .  of  course,  be  handle'll  by  ant  I .  li.atirig  its  effect  and 
calling  for  thrust  termination  when  V  ,  is  still  positive  by  the  required  amount. 
The  real  difficulty  comes  from  the  fact  that  there  is  some  uncertainty  attaw:hed  to 
the  residual  Impulse  due  to  valve  operating  times  and  other  factors.  A  typical 
number  here  might  be  151^  of  the  total  residuel  impulse,  or  4.8  fps.  Since  velo¬ 
city  uncertainties  of  this  site  are  not  tolerable  in  a  high  quality  guidance 
system,  liquid  engines  are  frequently  equipped  with  vernier  engines. 

Vernier  engines  are  small  liquid  or  solid  engines  which  are  capable  of  accelerating 
the  vehicle  at  0.1  -  0.2  g's.  When  a  vehicle  is  equipped  with  verniers,  the  main 
engine  is  shut  down  with  V  positive  by  an  amount  which  is  larger  than  any 
^  which  is  expected  from  a  residual  impulse.  The  vernier  lu  then  ignited  (If 
a  liquid  it  has  probably  been  on  with  the  main  engine)  and  its  thrust  is  terminated 
when  V  is  positive  by  an  amount  which  compensates  for  its  anticipated  residual 
impulse.  While  the  shut  down  procedure  is  the  same  for  the  small  engine,  the 
difference  Ilea  in  the  fact  that  the  velocity  uncertainty  is  reduced  by  the  ratio 
of  the  thrusts  of  the  main  emd  vernier  englnec.  For  the  0.1  g  vernier,  this 
means  that  the  velocity,  for  the  number  chosen,  would  be  reduced  by  a  factor  of 
100  to  O.OUfl  fps,  which  would  uoually  be  called  negligible.  For  medium  accuracy 
systems,  the  added  complexity  of  the  vernier  system  often  precludes  its  use. 


Pa^e  57 


Cutoff  Extrapolation  -  Even  when  a  vernier  ayateia  io  there  ntll] 

remains  the  problem  of  generating  the  thrust  termination  vith  sufficiently 

ranall  time  quantitation.  For  many  inertial  guidance  corapiiters,  position,  velocity 
and  veloclty-to-be -gained  are  computed  otce  every  one  half  r.ccond.  In  the  case  of 
a  0,1  g  vernier,  the  one  quarter  second  uncertainty  could  cause  a  0.8  fps  velocity 
error.  Shutlng  dovn  the  main  engine  with  a  time  uncertainty  this  large  wuuli 
also  require  the  use  of  a  vernier  with  longer  mean  burning  time.  In  order  to 
reduce  these  effects,  an  extrapolation  procedure  Is  used.  Consider  for  example 
a  system  where  V  is  generated  every  one  half  second  and  where  the  extrapola- 
tion  is  accomplished  every  one  sixteenth  of  a  second.  Then 


gx(N-*-n) 


V  V 

’'gxN  lb  V 


(117) 


where  V  „  Is  the  value  of  V  at  a  half  second  point  and  V  v  Is  the 

gxN  gx  ^  gx(N>n) 

value  n/lG  seconds  later.  Gufficlent  accuracy  is  usually  obtained  by  using  the 

standard  value  of  V  so  the  main  engine  Is  cutoff  when 

gx 


V  -  ♦  nK,  <,  K 
gxn  1  neco 


(118) 


and  the  vernier  engine  is  cutoff  when 


nK-  <  K  - — ' 
gxN  2  —  veco 


(119) 


where  K.  suid  K.,  are  standard  values  of  V  /l6  before  main  engine  and 
J.  <■  KXO 


.’cmier  cutoff  and  K 


^veco  computed  to  compensate  for  nominal 


values  of  residual  impulse  plus  on  uncertaiinty  in  the  case  oi'  the  main  engine. 

In  practice,  the  nominal  vernier  i>crlod  should  be  long  enough  to  damp  out  any 
steering  transients.  Usually  the  rco-  ired  length  is  in  the  20-30  second  region. 
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Systems  With  Solid  Hangings  -  Solid  propellant  enj^lnes  differ  from  liquid 
propellant  engines  in  thot  combustion  occurs  in  the  full  length  of  the  engine 
cane,  rather-THarf^ln  the  thrust  chamber  cnly.  Thmist  is  terminated  in  a  solid 
engine  by  bloving  holes  In  the  foi'vard  end  of  the  case.  This  has  the  effect 
of  reducing  the  thrust  to  z'*ro  in  a  very  short  period  of  time.  Even  vith  large 
engines,  the  residual  impulse  uncertainty  can  be  reduced  to  the  point  where 
verniers  are  not  required. 

Cutting  off  a  high  acceleration  engine  does  however  impose  some  rather  stringent 
requirements  on  the  cutoff  tine  granularity.  For  instance  for  a  ID  g  engine, 
cutting  off  to  -t-  1/32  second  would  result  in  a  velocity  uncertainty  of  10  fps. 

If  the  time  granularity  is  decreased  to  a  half  millisecond,  then  the  velocity 
uncertainty  is  decreased  to  0.16  fps.  Therefore  some  form  of  extrapolation  is 
required,  since  the  basic  guidance  computations  are  done  at  a  much  slower  rate. 
Here  it  turns  out  that  good  results  are  obtained  when  the  engine  is  cutoff  when 

V  =  V  ■*•/  a,„dx<K  (120) 

gx  gxo  /  -  meco 

^t 

o 

where  V  is  the  value  of  V  at  time  t  where  the  extranolation  mode  is 
gxo  gx  0 

begun  and  where  K _ _ is  a  constant  which  coTopensates  for  a  nominal  re  .si  dual 

meco 

Impulse.  The  extrapolation  assumes  that  and  the  integral  of  gravitational 

acceleration  change  a  negligible  amount  over  the  extrapolation  interval  which  is 
usually  true  for  intervals  of  the  order  of  a  half  second. 

Multiple  Thrust  Periods  -  As  has  been  pointed  out  previously,  more  complicated 

missions  such  as  placing  a  satellite  into  a  hlph  altitude  orbit  often  have  two 

or  more  thrust  periods.  D-urlnR  the  variou.s  thrust  period  ,  th^  ncminai  thrust 

direction  in  general  has  different  orientations  with  rerpect  to  V  ,  since 

gx 

the  computational  coordinate  system  usually  remains  fixed.  In  thla  ca.co  it  is 
necessary  to  resolve  either  V  or  the  steering  commands  themselves  from 

_  o 

Inertial  to  vehicle  coordinates.  This  can  be  done  either  in  the  computer  or 
by  means  of  an  analog  resolver  chain.  Cutoff  la  generated  on  the  basis  of  some 
major  component  of  V 


